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We study the pairing symmetry of a two orbital Ji — J2 model for FeAs layers in oxypnictides. We 
vary the doping and the value of Ji and J2 to compare all possible pairing symmetries in a mean-field 
calculation. We show that the mixture of an intra-orbital unconventional s^2y2 ~ cos(fca;) cos(fcy) 
pairing symmetry and a d^2_y2 ~ cos{kx) — cos{ky) pairing symmetry is favored in a large part 
of Ji — J2 phase diagram. A pure 53.2^2 pairing state is favored for J2 >> Ji. The signs of the 
dj,2_y2 order parameters in two different orbitals are opposite. While a small d^cy ~ sm{kx) sm{ky) 
inter-orbital pairing order coexists in the above phases, the intra-orbital dxy pairing symmetry is 
not favored even for large values of J2. 

PACS numbers: 



High temperature superconductivity (at 56K) has 
been recently reported in the rare-earth electron and 
hole-doped oxypnictide compoundsi'^'^'^'^'^. Prelimi- 
nary evidence^'®'^ suggests that the superconducting 
state in the electron-doped oxypnictides, like that in the 
Cuprates, has gapless nodal quasiparticle excitations and 
hence an unconventional pairing symmetry. A number of 
theoretical studies have predicted or conjectured differ- 
ent possible pairing symmetries, anywhere from p-wave 
to a TT-shifted s — wave^. 

A natural approach to the physics of oxypnictides is 
by drawing comparisons with that of the Cuprates. In 
the case of Cuprates, superconductivity is produced by 
doping a half-filled antiferromagnetic insulator. The an- 
tiferromagnetic exchange naturally provides for a singlet 
pairing amplitudeii, with equal mean-field critical tem- 
peratures for both d—wave and extended s — wave^^'^'^ in 
the ultra-Mott limit. Upon doping, due to the Fermi sur- 
face shape, the d-wave condensate has higher mean-field 
transition temperature than the extended s-wave one^^. 
The nearest neighbor hoping in the square lattice of Cu 
atoms is dominant and provides for a large Fermi surface 
with large effective mass in the Mott limit. 

The electronic properties of oxypnictides differ from 
those in Cuprates in several important ways. Most im- 
portantly, the undoped oxypnictides are metallic but 
their resistivity is strikingly high. They can hence be 
interpreted both as a bad metal or as a poor insulator, 
leaving open the question of whether a weak or strong 
coupling fixed point governs their physics. From the 
band structure point of view, barring the existence of 
un-physically strong crystal fields, it seems likely that all 
3d orbitals of the Fe atoms are involved in the low en- 
ergy electronic properties. Numerical results based on 
first principle calculations show the presence of small 
Fermi surfacesi^. In the unfolded Brillouin zone consist- 
ing of one Fe per unit cell, electron and hole pockets exist 
around the M and F, (tt, tt) points respectively. The mag- 
netic properties of the oxypnictides are also different from 
those of the Cuprates. Neutron experiments have shown 



that the magnetic structure in undoped oxypnictides is 
not a simple antiferromagnetic order— but rather a stripe 
spin-density wave with onset temperature of about 150K. 

The metallic behavior and the existence of Fermi pock- 
ets have led to proposals about the superconducting pair- 
ing symmetry and mechanism which originate in the 
weak coupling, itinerant limit^*^'^^. In this limit, triplet 
pairing is possible and has been predictedii. Numerical 
and analytic research suggests that the antiferromagnetic 
exchange coupling between Fe sites is stron g^^i^^i^" . Due 
to As-mediated hopping, an antiferromagnetic exchange 
coupling exists not only between the nearest neighbor 
(NN) Fe atom sites but also between next nearest neigh- 
bor (NNN) sites. The NNN coupling strength J2 is com- 
parable to the NN couphng strength Ji. The Ji — J2 
model provides for half-filled magnetic physics consistent 
with experimental neutron dataM. A nematic magnetic 
phase transition has been predicted in this model^i^, 
consistent with the experimental observation of a struc- 
tural transition preceding the spin density wave forma- 
tion. Therefore, the magnetic structure of the undoped 
oxypnictides suggests that the material is not far from 
the strongly coupling, Mott limit. 

In the present paper we obtain the superconducting 
mean-field phase diagram of a i — Ji — J2 model with 
the correct Fermi surface for the oxypnictide compounds. 
We predict that two kinds of intra-orbital pairing order 
parameters, an extended s-wave of the unconventional 
form Sx2y2 ~ cos cos ky or a dx2-y2 ~ cosfc^; — cos ky 
wave order parameter are the only possibilities in the 
Mott limit. A mixture of the intra-orbital unconventional 
Sx2y2 pairing symmetry and the dx2_y2 pairing symme- 
try is favored in a large part of Ji — J2 phase diagram. 
While a small d^y ~ sin^k^) sm(ky) inter-orbital pairing 
order coexists in the above phase, the intra-orbital dxy 
pairing symmetry is not favored even for large value of 
J2 in contradiction with the predictions of several pa- 
p3i- g20^22,23 ^j^g^^ j.gg^ analogy with the physics in 

cuprates. While dxy pairing would indeed be favored for 
J2 ^ Ji in the case of a large, single band, Fermi surface 
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FIG. 1: The phase diagram in Ji — J2 plane at 18% electron 
doping. 

(as in the cuprates)^, our calculation shows that, for the 
oxypnictides' Fermi surface, the dxy order parameter is 
the least favored. If the Fermi surface picture emerging 
from LDA calculations is correct, we can argue, on gen- 
eral grounds, that d^y pairing cannot compete with the 
extended s-wave we propose even if J2 is very large. If 
we consider a single band and treat a NNN J2 in mean- 
field decoupling, the superconducting transition temper- 
ature Tc is self-consistently determined by a Eliashberg 
equation 2Tc = J2Y.kU{k)f9{x{k,Tc)) where /(fc) is 
the pairing symmetry factor and g{x) — (with 

x{k,Tc) = '^*'2T~^ ) positive and peaked at the different 
Fermi surfaces. Hence Tc follows the maximum value 
of the pairing symmetry factor |/(A:)| close to the Fermi 
surfaces. Since in the unfolded Brillouin zone of oxypnic- 
tides the electron pockets are located around (0, tt), (tt, 0) 
and the hole pockets are located around (0, 0), (tt, tt), for 
the dxy symmetry pairing, the pairing symmetry factor, 
sin kx sin ky , is always small and dxy pairing symmetry is 
not favored even for large J2 . 

The model: We focus on a two-orbital per site model of 
the oxypnictides, with hybridization between the dxz and 
dyz orbitals. Although this description is only truly valid 
in the case of un-physically large crystal field splitting, 
we particularize to this model for analytic simplicity. We 
adopt the band structure proposed in Rcf.^^, which at 
first sight captures the essence the DFT results: 

T{k) = ( '^^^\~^' \ (1) 

where V^,^ = (4,,,fe,a' 4„,,fc,o-) is the creation operator 
for spin tr electrons in the two orbitals and the kinetic 
terms read: 

£x{k) — —2ti cos kx — 2^2 cos ky — At^ cos kx cos ky 
ey{k) = —2^2 cos kx — 2ti cos ky — 4^3 cos kx cos ky 

txy{k) — —At/^smkxSmky (2) 

The hoppings have roughly the same magnitude with 



h = -1.0,^2 = 1.3,^3 = -0.85, and U = -0.85. We 
find that the half-filled, two electrons per site configura- 
tion is achieved when /i = 1.54. The interaction Hamil- 
tonian contains three terms: the first two are an anti- 
ferromagnetic NN and NNN coupling between the spin 
of identical and opposite orbitals: 

-ff i = ^ Jiap{S{r, a) ■ S{r + 6i,P) - n{r, a)n{r + Si,l3)) 

(3) 

where S{r,a) — cl^ j,,^(j^j'Ca^r,'y' is the local spin opera- 
tor, n{r,a) is the local density operator, a,/3 are orbital 
index, i = 1,2, is the nearest neighbor and 62 is the 
next nearest neighbor. The third is a Hund's rule cou- 
pling of the spins on different orbitals, on the same site: 

H3 = -Y,JHS{r,a)S{r,a) (4) 

ra 

where a is the orbital complementary to a. This Hamil- 
tonian can be justified through different means. Hund's 
rule is known to play an important effect in Fe, but is 
usually neglected in the recent calculations on supercon- 
ductivity in oxypnictides. The antiferromagnetic Ji and 
J2 (both positive) are usually obtained from numerical 
calculations, although, in the Mott limit, they can be 
justified through a Hubbard-L/ Gutzwiller method. In 
the present case, numerical calculations predict a Hub- 
bard U for electrons on the same site and same orbital 
a factor of 2 larger than the Hubbard U' repulsion of 
electrons on the same site but different orbitals. In this 
case, the double occupancy constraint can be imposed 
by a product of the Gutzwiller projectors for the two or- 
bitals Pg = rijll - '^-^rfx.,T»^-^,dx.,i)(l " "-^d„,,T'^z,d^,a)- 
Simple perturbation theory leads us to Anderson ex- 

change 4(if/) V(C/+2 J^^)(^,,„5, ,^ — ni,c«ni,/3) where 
a, l3 are the two dxz and dyz orbitals. The hoppings in 
Eq[2] , give rise to NN exchange only between spins on the 
same orbitals, and NNN exchange between spins on both 
the same and opposite orbitals. However, other models 
for the spin-spin interactions, such as the one in Ref.'^-, 
which implicitly take into account Hund's rule by formu- 
lating the exchange in terms of spin-2 or spin-1 variables 
at each site, will also contain NN exchange between spins 
of opposite orbitals. Our mean-field solutions should be 
interpreted in the same spirit as the superconducting so- 
lutions of the original t — J model: at some value of the 
doping, the true undoped spin density wave groundstate 
disappears and gives way to a superconducting state^^. 

Mean-field solution in the absence of orbital crossing 
exchange: Keeping all of the above terms becomes ana- 
lytically intractable. We proceed with a two-step process: 
we first mean-field decouple the interaction Hamiltonian 
assuming that exchange takes place only between spins 
on the same orbitals. While physically incomplete, this 
model has the advantage of being analytically tractable, 
and exposes the un-competitiveness of the dxy order. We 
then numerically solve the full model in a superconduct- 
ing mean-field decoupling. The interaction term reads 
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'\'^a, — k,l'^'^; — k',l^ot.k','l 



with 



Vk,k' = -^J2±i'^os ± cos ky){cos k'^ ± cos k'y) (5) 



8J2 

N 



(cos kx cos ky cos k'^ cos k' + sin k^ sin ky sin k'^ sin k' 



with obvious pairing amphtudes in four channels ± 
2/2, xy and x^y^, Aa{k) = A^2+j,2 „(fc) + A3,2_j^2_„(fc) + 
Ax'^y'^,a{k) + A^y^aik), and 



cos fcx ± cos ky 



2Ji 

'1^ 



^(cosfc^ ± COS k'y)d{k') 



Ax^y^ (^) 

COS kx COS k,, 



Q T 



N 



sin kx sin A;,, 



where c?(A:') — {ca.^w ,iCa,k' ,])■ We use a — {1,2} to 



represent the orbital index [xz^yz). 

We decouple the interaction Hamiltonian with ex- 
change terms only between spins on the same orbitals 
in mean-field: Hm = Efc ^i^h 



A{k) 



/ ex(k) - 11 Ai(fc) f-xy{k) 

A|(fc) -e,(fc)-f/i -e.y(fc) 

ei:y(fc) ey{k) ~ ^1 A2(fc) 

V -exyik) A*{k) -ey{k)+fi. 



with ^'(fc) (ci,fc,T,cJ._fe i,C2,fc,T,4,-fc.i)- "^^"^ P^"^" 
ticularization to oxypnictides is present in the hop- 
ping terms, which couple different orbitals as in Eq[2]. 
A{k) can be diagonalized by an unitary transformatoin, 
U{ky A{k)U{k), and the Bogoliubov quasiparticle eigen- 
values El = —E2 and E3 = —E4 are given by 



(7) 



g2 

y 



Af + A2) ± Ji~el - ~el + A\- AlY + AeU{ix + ly? + (Ai - A^)^] (8) 



where tx,y = ^x.y — fJ-- The self-consistent gap and density 
equations are 



Ai(fc) = J2 yk.k'Ul^{k')Uim{k')F{Em{k')) (9) 

k' ,m 

Mk) = J2 yk,k'U:mik')U3mik')FiEm{k')) (10) 

k' ,m 

= 2^ C/(\,3)„(fc');7(i,3)m(fc')^(i?m(fc'))(ll) 



where F{E) is the Fermi-Dirac distribution function, 
F{E) = y+^s/et- To obtain the transition temperature, 
we linearize the self-consistent equation for small Ai, A2. 
After tedious algebra, we find the self-consistent equa- 
tions around Tc 



where 



W, = 



A2(fc) - ^kM'iWzik') - Wi{k')) 



{{ex-^lr-Ef)A, + elyAl 



(12) 



2\ex + ey-2fi\E,Jiely + iex-ey)^ 



:tanh(/3£;,/2) 
(13) 

with E, = E,{Ai = A2 = 0). 

The above equations can be solved numerically, vary- 
ing the doping fi and the value of Ji and J2. In FiglU 
we plot a phase diagram in Ji — J2 plane with 18% elec- 
tron doping. The phase on the left upper corner where 




FIG. 2: A three dimensional plot of the pairing weight Ws 
Wi as a function of (kx, ky) (electron doped) by setting Ai 
A2 = 1 in EqUl 



J2 > J2C 1.2 has pure extended s-wave pairing sym- 
metry Sx'^y^ phase. The phase on the right side, where 
Ji > Jic - 1.05, is a mixture of dx^-y^ and Sx'^+y^- 
The remaining large part of phase diagram is described 
by a phase with mixed Sx^y^ and dx2-y2 pairing order 
parameters. This part of phase diagram is believed to 
describe the real material since the estimated values of 
J2 ^ 0.5 and Ji ~ ij ^ 18,20 g^j-g sn^aUg]- than J(i,2)c- 
In this mixed state, the signs of the dx^-y^ order pa- 
rameters in the two orbitals are opposite. Namely, if 
Ai = acos(fc2;) cos(fcy) + b{cos{kx) — cos(fcy)), A2 = 
acos(fca;) cos(fcj,) — h(cos(kx) — cos{ky)). Moreover, we 
do not find a dxy solution in the entire parameter re- 



4 



gion. Time reversal broken states, such as s + id, are not 
favored either. 

The above results can be understood analytically. 
First, we can plot a pairing weight W3 — W\ as a func- 
tion of the Brillouin zone momentum (kx,ky) (Fig|5]) 
by taking A2 = Ai = 1 in EqlH The values of 
order parameters are determined by the pairing sym- 
metry factor function times this quantity. The dom- 
inant contribution is clearly around F, M and {tt,7t). 
The dxy order, in which the pairing symmetry factor, 
sinkx sinky, is peaked around (±7r/2, ±7r/2) has small 
overlap with the pairing weight and is not favored. Sec- 
ond, the mixing strength of two order parameters is 
determined by multiplying the two symmetry factors 
(/i; 12) of two order parameters and the paring weight: 
Efe/i(fc)/2(fc)(Vp3(fc) - Wiik)). It is easy to check, for 
a mixture of Sx2y2 and dxi-yi, i.e. /i = cos fc^; cos fcj^ , 
/2 = (cos kx cos/cy) that the summation has a large 
contribution from the Brillouin zone momentum around 
the electron pocket. The mixture strength of the other 
two order parameters ( Sx2y2 and Sx^+y^) is very small. 
This explains why the phase diagram is dominated by the 
mixture of 33,2^,2 and dx^^y-i . Finally, the difference of the 



relative sign between the Sx2y2 



and drr 



_y2 order param- 
eters in the two different orbitals is a result of the fact 
that exchanging kx to ky maps the xz to the yz orbital. 

The part of the phase diagram in Fig[T] with mixed 
Sx2y2 and dx2-y2 pairing becomes larger as the electron 
doping concentration is reduced: the mixing strength of 
Sx2y2 and dx2-y2 order parameters is (very slightly) in- 
creased due to the enhanced contribution around the M 
points. In Fig[3l we plot the transition temperature as 
a function of electron doping level at the fixed values of 
Ji — 0.25 and J2 = 0.5. On the electron-doped side, Tc 
is reduced by increasing the doping concentration. This 
is similar to Ref.^^ and it is, of course, around half filling, 
an artifact of the mean-field solution. The true ground 
state at half-filling is a spin-density wave^^ which gives 
way to a superconductor as the filling is increased. 

Solutions including orbital crossing exchange coupling 
and Hunds coupling: We now consider the orbital cross- 
ing exchange antiferromagnetic coupling, Ji;i2, J2;i2 and 
Hunds coupling Jh- In mean field, we can decouple 
this interaction in the particle-particle channel. The 
orbital crossing exchange coupling can be decoupled in 
four spin-singlet orbital crossing pairing order param- 
eters, A'(fc) = A',^2+y2{k) + A'^2_y2{k) + A'^2y2{k) + 
A'^yik). Hunds coupling can be decoupled to an on- 
site spin-triplet, orbital-singlet, order parameter. Ah = 

Efc(ci,fc,TC2,-/c,j, - C2,fc,T^i,-fc,i)- "^^^"^ ^^1^ 

Hamihonian is H'^ = j^k *(fc)^B(fc)*(A;) with B{k) = 
A{k) + SA{k) 
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FIG. 3: The pairing transition temperature as a function of 
the electron doping concentration at Ji = 0.25 and J2 = 
0.5. The dashed line indicates the region where the SDW 
competing phase takes over from the superconducting phase. 



We have an additional self-consistent equation: 
A'(fc) + A^ - ^ vl^,,u'^:,{k')uUk')F{EUk')) (14) 

k' ,m 

where the inter-orbital potential contains NN 



2Jia 
N 



'^j_{coskx ± cosfcy)(cosfc^ ± cosfc^). 



jj-^ (cos kx COS ky COS k'^ cos ky 



coupling 

a NNN coupling 
sin kx sin ky sin k'^ sin ky) and Hund's rule — The self 
consistent equations are solved numerically. We find that 
in the region where Jh Max{Ji, J2), Ah is extremely 
small. Hence, Hunds coupling has little effect on pair- 
ing symmetry. In the mixed S2;2j,2 and dx2-y2 phase, 
for Ji;i2 < Ji and J2;i2 % J2, the orbital crossing pair- 
ing order A' is zero within computing error except for 
dxy We find that a coexisting small inter-orbital par- 
ing order with dxy symmetry, A'(fc) = Aq sin(A;j;) sin(fcy). 
In Figdl we plot the result for the intra-orbital pair- 
ing order parameters 53,2^^2 and dx2_y2, and the inter- 
orbital pairing order parameter dxy as a function of 
J ~ Ji = J2 — Ji,i2 = •^2,12 when the chemical poten- 
tial is = 1.8 - corresponding to 18% electron doping. 
The result is a direct consequence of the dxy symmetry 
matching between the orbital-crossing pairing and the 
orbital-crossing hopping term. 

Discussion and Summary Preliminary experimental 
evidence suggests the presence of gapless nodal quasi- 
particle excitationa^i^i^ in oxypnictides. Our model pre- 
dicts a mixed intra-orbital order parameter Ai,2 — 
Ss{cos kx cos ky) ± Sd{coskx — cosky). In general, the 
mixed state is gapped, but with a very small gap. For a 
typical value of Sg = 0.2 and Sd = 0.1, the mixed state 
will have a gap of 0.025, which is around one fifth of the 
averaged superconducting gap in momentum space. The 
sign change of the d-wave order parameters between the 
two orbitals and the sign change of the s-wave order pa- 
rameter between the hole pockets and the electron pock- 
ets are interesting features of this state. The sign change 
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FIG. 4: The intra-orbital pairing order parameter, s^iyi, 
d^2_y2 and the inter-orbital pairing order parameter, d^y, as 
a function of J = Ji — J2 when the chemical potential is 
H = 1.8. 



can generate new physics, such as new bound states^! 



formed by impurity scattering. 

Although our prediction is based on a two-orbital 
model, we believe that the pairing symmetry predicted 
should be robust even if other orbitals are added. The 
pairing symmetry induced by the antiferromagnetic ex- 
change coupling is mainly determined by the structure 
of Fermi surfaces. As the Fermi surfaces in oxypnictides 
are located at F and M points, the dxy paring symme- 
try never wins over Sj;2j,2. Moreover, we find that the 
inter-orbital pairing is small even if the orbital crossing 
exchange is strong. 
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